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Recently Ciuchini, Pierini and Silvestrini proposed a method for constraining CKM parameters 
in B — > Ktttt and B s —* Ktttt through phase measurements of amplitudes involving 7 = 3/2 K*tt 
final states. We show that complementary information on CKM parameters may be obtained by 
studying the phases of AJ = 1 B — » (K*ir)i = i/2, B s — > (K* K)i=i and B s — > (K*K)i=i amplitudes. 
Hadronic uncertainties in these constraints from electroweak penguin operators Og and Oio, studied 
using flavor SU(3), are shown to be very small in B — > Knn and B s — > Ktttt and somewhat larger in 
B s — » KK-n. The first processes imply a precise linear relation between p and r], with a measurable 
slope and an intercept at fj — involving a theoretical error of 0.03. The decays B a — * A"7T7t permit a 
measurement of 7 involving a theoretical error below a degree. We note that while time-dependence 
is required when studying B° decays at the T(4S I ), it may not be needed when studying B s decays 
at hadronic colliders. 



I. INTRODUCTION 

Recently a method has been proposed by Ciuchini, 
Pierini and Silvestrini P, for determining Cabibbo- 
Kobayashi-Maskawa (CKM) parameters in three body 
B — > Ktttt and B s — » _ftT7T7r decays. The proposed 
method is reminiscent of early suggestions for determin- 
ing 7 using rates and asymmetries in two body decays 
B — > Ktt §| E IE El and B s — > Ktt 7} . A unique fea- 
ture of the new method is being able to measure through 
interference in the Dalitz plot relative phases between 
quasi two-body decay amplitudes for Br s \ — > K*tt and 



B 



K*tt. This is similar to a proposal for measuring 



relative phases among B — > p7r amplitudes by studying 
the Dalitz plot for 5° -> tt+tt^tt 0. When neglect- 
ing electroweak penguin (EWP) contributions, the rela- 
tive phase between a combination of decay amplitudes 
describing Br s \ — > (K*tt)j =3 /2 and a corresponding com- 
bination of -B( s ) amplitudes determines the weak phase 7. 
A small hadronic uncertainty caused by EWP amplitudes 
was estimated, based on factorization and assuming cer- 
tain input values for B-to-light-mesons form factors pj. 

In the present paper we propose extending the method 
to AJ = 1, 1(K*ir) = 1/2 amplitudes in the above decays 
and to I = 1 amplitudes in B s — > K*K and B s — > K*K. 
We use flavor SU(3) to study theoretical uncertainties 
caused by EWP contributions, suggesting a way for re- 
ducing these uncertainties. The resulting theoretical pre- 



cision in determining CKM parameters in B — > Kttt: and 
B s — » Ktttt is shown to be very high, essentially at the 
level of isospin breaking corrections. This happens be- 
cause the method is based primarily on isospin symmetry 
considerations, while flavor SU(3) is used only to estimate 
uncertainties from a subset of small EWP contributions. 

In Section II we analyze B — ► K*tt, B s — * K*tt and 
B s — * K* K(K* K) decays in terms of isospin ampli- 
tudes, selecting several ratios of Br s \ and B/ s \ isospin 
amplitudes which can be used to determine 7 in the ab- 
sence of EWP contributions. Section III studies the ef- 
fects of EWP amplitudes, turning the determination of 
7 into a generic constraint on CKM parameters. The 
constraint involves an uncertainty from a ratio of two 
hadronic matrix elements of (V — A) current-current op- 
erators. Flavor SU(3) calculations show that this ratio 
is small for judiciously chosen combinations of isospin 
amplitudes in B — > K*tt, vanishes approximately in 
B s —>■ K*tt in the isospin symmetry limit, and is larger 
in B s — * K* K(K* K). This implies precise constraints 
on CKM parameters from knowledge of amplitudes and 
their relative phases for B — > K*tt and B s — * K*tt. 

Section IV discusses measurements of these quasi two- 
body decay amplitudes and of B — + K*K (K*K) in three 
classes of three body decays, B — > Ktttt, B s — * Ktttt and 
B s — > KKtt decays, respectively. We point out that 
measuring a relative phase between the amplitudes for 
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K*+tt- and B° -> K*-tt+ in B° -> K s tt+tt- 
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produced in e + e~ collisions at the T(AS) requires time- 
dependence. In contrast, no time-dependence may be 
needed for a similar measurement in B s —> Kstt + tt~ per- 
formed at hadronic colliders if a width difference in the 
B s system is measured. In order to obtain a most precise 
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determination of CKM parameters, we propose applying 
amplitude analyses to the entire B — *■ Kmr class, using 
isospin amplitudes as variables. Section V concludes with 
several remarks about the implementation of this method 
and its sensitivity to physics beyond the Standard Model, 
comparing it with two other methods for determining 7 
in B and B s decays. 



II. ISOSPIN DECOMPOSITIONS AND 7 
WITHOUT ELECTROWEAK PENGUIN TERMS 

The cleanest method for extracting the weak phase a 
or 7 in AS — and AS = 1 charmless hadronic B decays 
stems from applying isospin symmetry to these decays, 
eliminating the effect of QCD penguin amplitudes which 
transform in these processes as AI = 1/2 and AI = 0, 
respectively ;9| . We will now discuss separately the four 
cases, B -> K*ir, B s -> K*ir, B s -> K*K and B s -> 
K*K in terms of their isospin amplitudes. 



A. B — > K*ty 

In strangeness changing decays of the type B — > Ktt 
or B — > X*7r, where K* denotes any kaon resonance, 
if* (892), K*(U30), iq(1430), . . . , the four physical am- 
plitudes for B° and B + are decomposed into three 
isospin-invariant amplitudes [3, 



-A(K*+tt-) = B 1/2 - A 1/2 - A 3/2 , 
A(K*\+)=B 1/2 + A 1/2 + A 3/2 , 
-V2A(K*+tt°) = B 1/2 + A 1/2 - 2A 3/2 
V2A(K*°ir°) = B 1/2 - A 1/2 + 2A 3/2 



(1) 



Here we use a phase convention in which a minus 
sign is associated with a u quark in a meson. The am- 
plitudes B and A correspond to AI = and AI = 1 
parts of 7i ff, respectively. Their subscripts denote the 
isospin of the final K*ir state. Here and elsewhere we will 
denote by B isospin amplitudes obtaining contributions 
from QCD penguin operators, and by A other ampli- 
tudes. Our study will focus on the latter. 
The amplitude quadrangle relation, 



3A 3/2 = A(K*+tt-) + V2A(K*°ir°) 
= A(K*°tt+) + V2A(K* + ir ) 



(2) 



defines A 3 / 2 as one diagonal of the quadrangle, while 
A\i 2 is given by 

6A 1/2 = A(K*+tt-) + 3A(K* a TT+) - 2V2A(K*°ir°) 
= 3A{K* + n-) + A(K* a 7T+) - 2V2A(K*+tt°) . 

(3) 

The two AI = 1 amplitudes, A 3 / 2 and A 1 / 2 , do not 
contain a QCD penguin contribution and would carry 



a single weak phase 7 if EWP contributions could be 
neglected. Here we will proceed under this assumption, 
postponing a discussion of the effects of EWP amplitudes 
to the next section. Denoting amplitudes for charge- 
conjugate initial and final states by A, and defining two 
ratios of amplitudes, 



Ri 



A! 



I = 1/2, 3/2 , 



the phase 7 is determined by 
1 



$/ = --arg(-Rj) = 7 



(4) 



(5) 



Note that although the ratios Rj do not depend on the 
magnitudes of Aj in the limit of vanishing EWP contri- 
butions, an extraction of 7 requires measuring both the 
magnitudes and the relative phases of physical B — > K*n 
amplitudes and their charge-conjugate. 

The ratio R 3 / 2 was studied in [l| (where it was denoted 
by R° — R ± ) while the ratio R\/ 2 studied here provides 
independent information on CKM parameters. 



B. B„ 



K*TT 



The isospin decomposition of the two AS = B s 
K*tt decay amplitudes is: 



A s {K*+ir-) = A s 3/2 -V2B s 1 
A s {K*°ir°) = V2A» 3/2 + B 



1/2 



(6) 



'1/2 



where the superscript s denotes B s instead of B Q and 
subscripts denote the isospin of both the transition oper- 
ator and the final K*ir state. Since in AS = decays the 
QCD penguin operator behaves as AI = 1/2 it is con- 
tained only in -Bjy 2 - On the other hand, the amplitude 

3A' /2 = A s {K*+ir-) + V2A s (K*°tt°) (7) 

is pure tree when neglecting EWP contributions, thus 
providing information on 7. Defining a ratio of B s and 
B s amplitudes (denoted Rd in |2j), 



I? 



.4 



3/2-^ 



3/2 



3/2 



one now has 



3/2 = -,arg(i^ /2 ) =7 



(8) 



(9) 



B a — » K*K and B s 



K*K 



These AS = 1 decays involve two independent pairs of 
isospin amplitudes, 



a s (k*+k-) = a s 1 + b s , 

A s (K* Q K )=Al-B° , 



(10) 
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and 

A S (K*-K+)=A' 1 S + B' S , 
A s (K*°K a ) = A'i - B' S . 

Thus, one has 

2A\ = A S (K*+K-) + A S (K*°K°) , 
2A'l = A S (K*-K+) + A S (K*°K°) . 



(11) 



(12) 



Defining for each of these processes a ratio of B s and B s 
amplitudes, 



AS A's 
ps _ A l p/s _ A l 

1 ~ If ' Kl - A'f ' 
one obtains two new independent equations for 7, 

$? = -iarg( J Rf)= 7 , 

$i s = -^arg«) = 7 . 



(13) 



(14) 



III. CKM CONSTRAINTS INCLUDING 
ELECTROWEAK PENGUIN AMPLITUDES 

In section ITT1 we have neglected AS = 1, AI — 1 and 
AS = 0,AJ = 3/2 EWP contributions. In this limit 

measurements of the ratios Rj in JSJ, 10 and i|14|) 
determine 7. The inclusion of EWP operators modi- 
fies these relations since these operators involve different 
weak phase than the tree operators. These effects are 
important in the AS* = 1 relations JSJ and l|14|) . where 
EWP contributions are CKM-enhanced, and are negligi- 
ble in the AS = relation JUJ). We will first obtain a 
general constraint in the (p, fj) plane pd| following from 
fixed values of <f>\ s) = -i a rg(i?^ ) ). 

Let us study the effect of EWP operators on obtain- 
ing CKM constraints in AS = 1 decays. The dominant 
(V-A) EWP operators, O s 9 and OJ , in the AS = 1 effec- 
tive Hamiltonian \l2^ are related to current-current oper- 
ators, Of = [s6]v_yi[uM]y_A and = [ub]v-A[su]v-A, 
through operator relations, 



^9,10 = 2^1-2 + [°P erators witri = 0] 



(15) 



Neglecting EWP operators, O7 and Of, involving small 
Wilson coefficients, the AI — 1 part of the AS* = 1 weak 
Hamiltonian can be rewritten as 



H 



AI=1 



I K C- 



-X!C 



:A?C; 



EWP \ Q AI=1 



(16) 



where A s 



AI=l 



M t) = V u(t)b V u(t)s, O^ 1 ^ U°l ± °2), and 
C± = Ci±C2, C^ WP = CgiCio are sums and differences 
of Wilson coefficients. 



Terms in (I16|) involving Cj? WP introduce in AI = 1 
amplitudes a weak phase different from 7, with coeffi- 
cients depending on hadronic matrix elements for O ai=1 
and 0+ I=1 . Using a relation between Wilson coefficients 
which holds up to 1% corrections |l2(, 



£<EWP (jEWP 



C 4 



o_ 



(17) 



one obtains a generic expression for the four ratios 
Ri/2,Rz/i,Rt and R[ s in Eqs. @ and (|T3J), 



R (s) _ c -2if 7 +arg(l + >c)1 1 + C l r I 



Here we define 



1 - K 
1 + K 



1 + c K r>' 



3 Of WP Af 
2 + A?, 



(Mo_o^ =1 |i? 



(18) 



(19) 



(20) 



The parameter k depends only on calculable Wilson 
coefficients and on CKM parameters. In order to il- 
lustrate the sizable shift in <E>^' = — |arg(i?^) caused 
by this parameter alone, we use the central values for 
CKM parameters ^3 an d next to leading order (NLO) 
values for Wilson coefficients at p = m b = 4.8 GeV, 
d(m b ) = -0.178, C 2 (m b ) = 1.079, C 9 (m b ) = -0.0102, 
O10 = 0.0017. We find 



X s 

K = -1(1.404 ± 0.038) x 10~ 2 



K 
= -0.35 



(21) 



0.56i 



where the error in the brackets corresponds to varying 
the scale p in the NLO Wilson coefficients in the range 
m b /2 < p < nib- In the second line we give the result for 
central values of Wilson coefficients and CKM elements. 
This value of k translates into arg(l + k) = 41°. 

(s) 

A nonzero value of the parameter r\ leads to an ad- 
ditional shift in ^ { j s \ given by -±arg[(l + c*r| s) )/(l + 

c K r\ s ^}. A given value of the observable can be shown 
to imply the following constraint in the {p, fj) plane (we 
use A = 0.227): 



V+(p + C)t 

(p + c)-fjt 



= tan $ 



Here we define 



_ 3 C™ p 1 _ A 2 /2 



Oee - 



2 4 



1 + M- 



A 2 



= -0.27 



t = 



v /(l + t 2 + )(l+< 2 _) 



(22) 



(23) 



(24) 
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FIG. 1: la constraint in the p — fj plane (two almost vertical 
parallel black lines) from a precisely measured $3/2 in B — > 
Kith, taking values for r 3 / 2 as in (l.'K it . All other constraints 
are taken from |l3j|. 



complex phase is related to 7, is of order a few percent 
[see Eqs. (JEJ) and l(2T)l]. 



4\ 



2C+ \\t\ 



= (1.40 ±0.04) x 10 



(27) 



-2 Sm 7 

sin/3 
k') on 7 is very 



Thus, the dependence of the shift arg(l - 
small and calculable in terms of 7. 

Another case where r = holds in a symmetry 
limit is i? — !■ (Kn)j =3 /2, where the K and it mesons 
are in an S-wave and must be in a symmetric SU(3) 
state 

HQ- This SU(3) argument does not hold in B 



and B s 

» 



1/2,3/2 



(K*K) 



1=1 



nor does it hold in B s — ► {K*K)i=\ 
We will study now the values of 



for these decays within flavor SU(3). Theoretical er- 
rors in these values lead to uncertainties in the resulting 
CKM constraints. 



Ratios ri and CKM constraints in B — > K*n 



(p 2 + fj 2 ~ C 2 )Im(r[ s) ) T 2Cf]Re(rY > ) 

(p + C) 2 + Tf + {p 2 +fj 2 - C 2 )Re(4 s) ) ± 2Cfjlm(r { j S) ) ' 

(25) 

For a small value of and for a given value of the 

(s) 

observable <&} , one obtains the following constraint, 



(s) 

T) = tan $} ' 



C 



l-2Re(rj s) ) 



+ 0(r[ s)2 ) . (26) 



This describes a straight line in the (p, rj) plane (cf. 
Fig- 01, with a slope tan^^ and an intercept po = 
-C[l - 2Re(r[ s) )] at fj = 0. A theoretical error <5r|' S ' ) in 



^ translates into an uncertainty of ±2|C|<5r} f ''' in the hi- 
ts) 

tercept po but no uncertainty in the slope tan <&} which 

(s) 

is measured through the ratio i?} . Assuming for illus- 

(s) 

tration a negligible value of r\ , one may estimate the 
slope required in the Standard Model by choosing central 
values of (p,fj) from a CKM fit [T^, p = 0.20, fj = 0.34. 

(s) 

This implies a slope tan<I>j = —5.0, which is quite sen- 
sitive to the value of . 

A similar treatment of EWP contributions can be ap- 
plied to AS = processes involving CKM factors ^uy 
In the isospin symmetry limit the AS = part of 0_ is 



(«) 



pure AI = 1/2, hence O 



AJ=3/2 



0. Consequently, in 



B s — > (^*7r)/=3/2 one has r|, 2 = 0. 1 This implies that 
there is no hadronic uncertainty in the CKM constraint 
from the ratio Ry 2 m B s — * K*w, aside from tiny cor- 
rections from the operators O7 and which we have 
neglected. The parameter k' in AS = decays, whose 



1 This observation has been overlooked in Ref. 0- 



The two ratios Ri in Eq. (JIJ, providing independent 
pieces of information on 7, are given by Eq. I|18[) with 
r = r 1 given by 



ri 



((K*7r) I \C-Q A1 = 1 \B) 
{{K**)i\C+0?=i\B) 



1/2,3/2. (28) 



The ratio r 3 / 2 was estimated in [lj, based on factoriza- 
tion and assuming certain input values for _B-to-light- 
mesons form factors. Here we wish to present a differ- 
ent approach based on flavor SU(3) for calculating both 
r 3/2 an d r i/2- Using flavor SU(3) rj may be calculated 
from tree-dominated strangeness-conserving B decays, 
which arc CKM-enhanced relative to tree amplitudes in 
B — ► K*tt, and which have already been measured. Fur- 
thermore, one may apply Eqs. I|18|l and l|28|l to \ (K*tt)x), 
an arbitrary superposition of / = 1/2 and 7 = 3/2 K*n 
states. The corresponding ratio of hadronic matrix el- 
ements will be denoted rx- One is searching for a lin- 
ear superposition of isospin states which leads to a small 
value of rx in order to obtain a small uncertainty in CKM 
parameters. 

The operators 0+ I=1 and 0^ I=1 transform as 15 and 
6 representations of SU(3), respectively, while a gen- 
eral K*tt state is a combination of 8s, 8a, 10, 10 and 
27 O E3. Thus, the numerator in rx involves in 
general a linear combination of three reduced matrix 
elements, (8 S |6|3), (8 A |6|3), (10|6|3), and the denomi- 
nator involves a combination of four matrix elements, 
(8 S |15|3), (8 A |15|3), (T0|15|3), (27|15|3). The same re- 
duced matrix elements occur in AS = amplitudes. One 
is seeking two sums of AS* = amplitudes which are 
given by the same two combinations in the numerator 
and denominator of rx ■ 

The case of r 3 / 2 is particularly simple since the state 
{K*Tt)i=3/2 contains only two pieces transforming as 10 
and 27 of SU(3). The transformation properties of the 
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TABLE I: Branching ratios for B — > pn and B — > K*K decays, in 
units of 10 — 6 , taken from Ref. |22| unless quoted otherwise. 





IvlULlt; 


JDI dlldllllfJ, IXclLlL) 


B+ 




°- ' — i.i 




P+7Y° 






K^K+ 


< 5.3 


5° -» 


p ± n^ 


24.0 ±2.5 




p + TV~ 


16.0 ± 2.0 a 




p~1T + 


7.6 ± 1.3 6 




J> o 

P 7T 


i qo+0.56 
1 -°' i -0.55 






< 1.9 



"We take an average of 19.5 ± 5. [p and 15.3 ± 2.2 El. 
6 We take an average of 9.6 ± 3.4 l23l and 7.3 ± 1.4 |24I. 



operators imply that only the 10 and 27 pieces con- 
tribute to the numerator and denominator, respectively. 
Thus, r 3 / 2 is proportional to a ratio of corresponding re- 
duced matrix elements, (10|6|3)/(27|15|3). The numeri- 
cal coefficient multiplying this ratio can be read off SU(3) 
Clebsch-Gordan tables in (after translating into our 
phase convention). These tables can also be used to ex- 
press (10|6|3) and (27|T5|3) in terms of AS" = ampli- 
tudes, 

r 3 /2 = 

[A{p+n°) - A(p°n+)] - V2[A(K*+K°) - A(K*°K+)} 



color factor, 



A(p+TT°) + A(p°n^ 



(29) 



This expression can be simplified by neglecting the 
AS = QCD penguin amplitude given by the second 
term in the numerator, and by assuming that the strong 
phase difference between the two amplitudes in the re- 
maining term is small, as this phase is expected to be 
suppressed by l/rrib and a s (rrib) 0,0, H3- This is sup- 
ported by studies of QCD penguin amplitudes (including 
charming penguins) in B ~ * pn which have been found to 
be small, with a penguin-to-tree ratio of about 0.2 |2lj . 
Signs of color-allowed amplitudes are assumed to be given 
by factorization. Using branching ratios given in Table 
I, one finds 



r 3 /2 



_ \y/B(p+n°)-y/B(J>n+)\ 
^B{p+n*) + ^B{p*n+) 
= 0.054 ±0.045 ±0.023 . 



(30) 



The first error is caused by experimental errors in B — > 
pn branching ratios. The second error, due to SU(3) 
breaking, is calculated by allowing an uncertainty of 30% 
in each of the reduced matrix elements entering the phys- 
ical amplitudes. 

The value Ip^Hjl . obtained by applying SU(3) to B — > 
pn branching ratios, may be compared with an estimate 
based on naive factorization 1] in which we include a 



\ BK" 



C- (1 - l/N c ) (UAfi 
7-3/2 C + (1 + 1/N C ) (f n A$ K * 
= 0.012 ±0.083 . 



+ f K *F Q B ") (31) 



We used the following values for decay constants and 
form factors = 131 MeV, f K * = 218 ± 

4 MeV, F B7r = 0.28 ± 0.05, Af K * = 0.45 ± 0.07. Note 
that naive factorization may be a reasonable approxima- 
tion because the ratio r 3 / 2 defined in l|28(l does not involve 
QCD penguin contributions. 

A bound on the error in r 3 / 2 caused by neglecting a dif- 
ference of two AS = QCD penguin amplitudes in 129|) 
can be obtained in terms of upper bounds on branching 
fractions for B+ -> K*+K° and B+ -> K*°K+. Al- 
though current upper bounds are not very useful (see see 
Table I), we expect the bounds to improve in the future, 
such that the error caused by neglecting these terms will 
be at most at the level of SU(3) breaking. 

The error in r 3 / 2 affects the CKM constraint <|26[) 
through the term involving Re(r 3 / 2 ). The error from 
neglecting a strong phase difference between Ap+^o and 
Ap0 7r + is expected to be very small, since Re(r 3 / 2 ) de- 
pends quadratically on this phase. This is gratifying since 
the size of 1/rrib sup pressed strong phases cannot be reli- 
ably calculated |18l l25j . Information on the above phase 
is provided by the isospin pentagon relation 



A{p + n°) + A{p°n + ) = 

±=(A(p + n-) + A{p-n + )) + V2A(p°n°) 



(32) 



Relative phases between amplitudes on the right-hand- 
side can be measured through a Dalitz plot analysis 
of B° — > n+n~n° 0, |2jj. Assuming that phases be- 
tween B — ► pn amplitudes can be neglected, and u sing 
branching ratios from Table I and a lifetime ratio |22| . 
r+ /r = 1.076 ± 0.008, Eq. reads 6.01 ± 0.27 = 

6.69 ± 0.38. This agreement shows that relative phases 
between B — > pn amplitudes are not large. Assuming 
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in contrast a negative sign for the color-suppressed am- 
plitude A(p°ir°), for which factorization does not hold, 
would imply 6.01 ± 0.27 = 2.86 ± 0.38 which is badly 
violated. 

The value of r 3 / 2 in l|30l) can now be substituted in 
Eq. H26|) ■ The resulting linear constraint in p — fj plane 
is shown in Fig. ^ assuming a precisely measured value 
for $3/2 ■ The current error in r 3 / 2 translates into a very 
small error of in the intercept where fj = 0, po = 0.24 ± 
0.03, but no theoretical error in the slope which is given 
by a value measured for tan$ 3 / 2 . The small error in the 
intercept, partly from SU(3) breaking in r 3 / 2 , is linear in 
the uncertainty in r 3 / 2 , and may be reduced only slightly 
by measuring more precisely B — > pit branching ratios. 

The calculation of T\ii proceeds in a similar manner 
to the calculation of r 3 / 2 , leading to a larger value of 
order one. Instead, one may search for a superposition of 
1 = 1/2 and 7 = 3/2 K*ir states for which rx is small. 
Using 



\(K*n) x ) = -j={\I = 1/2) - 2\I = 3/2)) , (33) 



we find 



rx 



A(p+ir°) - 2A{p°tt+) + V2A{p°tt°) + N X (K*K) 

A{p+n°) + V2A(p-ir+) + V2A(p n°) + D X (K*K) ' 

(34) 

AS* = QCD penguin and annihilation amplitudes in 
the numerator and denominator. 



Nx{K*K) =[2A(K*°K+) - 3A(K*+K°) 

- A(K*°K°) + A(K* + K~)] /y/2 

D X {K*K) =[A(K*°K Q ) - A(K*+K°) 
+ A{K*+K-)]/V2 , 



(35) 



are expected to be no larger than SU(3) breaking correc- 
tions and will be neglected. Assuming small strong phase 
differences between B — > pix amplitudes, and using mea- 
sured branching ratios in Table I, we find, 



rx 



-0.068 ± 0.057 ±0.044 



(36) 



The first error originates in experimental errors in B — » 
pir branching ratios, while the second error is calculated 
assuming 30% SU(3) breaking in reduced matrix ele- 
ments. The central value of rx and its error are com- 
parable to r 3 / 2 . As in the latter case, this translates to 
a very small error in the intercept, but no error in the 
slope of the linear relation (|26|l between p and fj provided 
by a measurement of $>x- We define Rx and <I>x as in 
Eqs. |@} and 10 using the K*ir state defined in (O, 



R > 



/2 



2A 



3/2 



2A 



3/2 



Ai/2 

<S>x =--a,rg(R x ) 



(37) 



The result ^ may be compared with an estimate 
based on naive factorization, 



rx = 



C- (1 - 1/A e ) (2f n A* K ' - f K -F B «) 

C+ (1 + 1/JV C ) QUA*?*' + f K *F B ") (38) 

-0.22 ±0.07 , 



where the error reflects only errors on the assumed form 
factors. Our result l|36|l using flavor SU(3) agrees within 
uncertainties with this more crude approximation which 
gives a somewhat larger central value. 



B. Ratio r| ; 2 and determining 7 in B s — > K*n 

As we have shown using isospin symmetry alone, the 
parameter r in B a — > K*it vanishes, r|/ 2 — 0> because 
the AS" = part of <3_ is pure AI = 1/2. The small 
parameter k! introduces a small shift arg(l + k') in 
away from 7. Since the shift is calculable in terms of 
7 [see Eq. I|27p]. the theoretical error in determining 7 
using these processes is below one degree. 

Note that measuring 7 in these processes, and using 
B — > Kttt: for constraining the point (p, fj) to lie on a 
straight line with measured slope and intercept, fixes the 
apex of the unitarity triangle as the point where the two 
straight lines intersect. Thus, in principle, B — > Kmr and 
B s — * Ktttt alone determine the shape of the unitarity 
triangle. 



C. Ratios rf,r'i and CKM constraints in 

B s -» K*K, K*K 

In the presence of EWP contributions the two ratios 
i?f and R'i defined in (|13H are given by Eq. JTBJ) with 



(K*K\C-0 AI=1 \B s 
(K*K\C + 0^\B s 
(K*K\C-0 AI=1 \B s 
(K*K\C+0 AI = l \B s 



(39) 



We use SU(3) tables in Ref. [Tjj to express these ratios 
in terms of AS = decay amplitudes for nonstrange B 
mesons, 

s _ A{p+tt-) + A{p-Ti+) - V2A(p°t:+) + Nx{K*K) 

Tl ~ A(p+n-) - A(p-n+) + V2A(p a Tr+) + Dx{K*K) ' 

la _ A(p-TT+) + A(p+tt-) - V2A(p+Tr a ) + N[{K*K) 

Tl _ A(p~tt + ) — A(p + n~) + \/2A(p + tt®) + D[(K*K) ' 

(40) 

Penguin and annihilation terms in the numerators and 
denominators, 

Ni(Dt) = ±A(K* Q K+) t A{K*-K+) + A(K*°K°) , 

N[(D[) = ±A(K*+K°) t A(K* + K~) + A(K*°K°) , 

(41) 
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will be assumed to be smaller than SU(3) breaking cor- 
rections. 

Disregarding phase differences between B — > pit ampli- 
tudes which have a very small effect on rf and r'* [as we 
argued in obtaining (|30[) ]. using measured branching ra- 
tios in Table I, and estimating errors from SU(3) breaking 
as explained above, we have 



0.52 ±0.10 ±0.18 
0.70 ±0.21 ±0.41 



(42) 



Comparing these values with an estimate based on naive 
factorization, we find agreement again, 



C_ (1 - 1/N C ) 
'C+ (1 + 1/JVc) 



0.70 



(43) 



We do not expect the errors in (|42|l to improve by 
reducing errors in B(B — > pir), as SU(3) breaking intro- 
duces a comparable uncertainty. The values of rf and 
r'i can be substituted in Eqs. (|22II - H25|I to obtain con- 
straints in the (p, fj) plane, for measured values of $f 
and <$>'*. The larger errors in (|42|l in comparison with 
those in (|30|l and i|36|l imply larger uncertainties in these 
constraints than in those following from <t> 3 / 2 and §x ■ 



IV. MEASURING MAGNITUDES AND PHASES 
FOR QUASI TWO-BODY DECAY AMPLITUDES 

As shown in the previous two sections, new con- 
straints in the (p, 77) plane can be obtained within each 
of the three classes of quasi two-body decay processes, 
B -> K*n, B s -> K*n, and B s -► K*K, K*K and 
their charge-conjugates. This requires measuring both 
the magnitudes of the amplitudes in a given class and 
their relative phases. This can be achieved through am- 
plitude analyses of charmless three-body decays which 
we discuss now. 

A three-body B (or B s ) decay amplitude into a fi- 
nal state /, which is a function of two Dalitz variables, 
Si2j S13, is expressed as a sum of several Breit-Wigner res- 
onant contributions and a non-resonant term. Resonant 
contributions are given by complex constant amplitudes 
Aj multiplying Breit-Wigner functions (S12, S13), 

while the nonresonant amplitude Anr may vary in the 
S12, S13 plane, 



BW 



(*i2,*is) ■ (44) 



The corresponding amplitude A(s 12, S23), for three-body 
B (or B s ) decays into a charge-conjugate state /, is given 
in terms of an amplitude .Anr and a set Ai correspond- 
ing to charge-conjugate resonances. In general, one has 
^nr 7^ ^NRj Ai ^ Ai as each amplitude may involve two 
weak phases and two different strong phases. Direct CP 
violation in a particular resonant or non-resonant channel 
would be implied by \Ai\ ^ \A t \ or |^4 N r| ^ I^nrI- 



Fitting the event distribution for three body B (or B s ) 
decays to the squared amplitude (|44|l permits determin- 
ing the magnitudes of A4 and their relative phases. A 
relative phase between two resonant amplitudes is di- 
rectly measurable when the two resonances overlap in the 
Dalitz plot. This relative phase can also be measured 
when there is no overlap between the two resonances, 
but each of the two resonances overlaps with a third res- 
onance. Alternatively, a phase between two resonance 
amplitudes can be measured through their interference 
with the nonresonant amplitude Anr- 

We will be interested primarily in relative phases be- 
tween amplitudes associated with K meson resonant 
states. Charmless three-body decays involving 7r + 7r~ or 
K + K~ obtain also contributions from cc resonant states, 
which involve relatively small rates and are expected to 
lead to sizable CP asymmetries [13, IH, ■ 



A. B — > Knn 

We start this discussion with the decays B — > Kirir 
which are currently the most feasible ones among the 
three classes studied in this paper. Amplitude analyses 
of B — > Ktttt, for both charged and neutral B mesons, 
have been performed by the Belle and Babar collabora- 
tions. Decays B + — > K + ir + ir~ have been studied by 
both Belle [3(j and Babar [3l]]. An amplitude analysis 
was made by Babar H2 for B° -> K+tt-tt° jH, and by 
Belle 34] for B° — > KsTT + n~ . The first two processes are 
self-tagging whereas the third decay involves final state 
which is not flavor specific. These measurements have 
already provided some useful information which is rele- 
vant to our proposed study. We note that these studies 
have averaged over the above processes and their CP- 
conjugates. The proposed study requires separate ampli- 
tude analyses for B and B decays. 

The process B + — ► K + tt + tt~ gave information about 
the magnitudes of amplitudes for B + — > K*° (892)n+ and 
B+ -> K*°(1430)7r+ and their relative phase [M Hl|. 
The statistical error in the measured relative phase is at 
a level of 10° which is encouraging. However, this three- 
body decay provides no information on a relative phase 
between two B — ► K*(892)ir amplitudes where pairs of 
K* and tt have different charges. 

The decay B° — > K + iy~tt° is more interesting in our 
context, since it measures the magnitudes of A(B° — > 
K*+n-) and A(B° -> X*°7r°), for both K*(892) and 
.Kg (1430), as well as the three relative phases among 
these amplitudes. Errors in the measured phases are at a 
level of 40° j32 • It would be useful to understand the ori- 
gin of this large error in order to reduce it in future stud- 
ies of this process, and to perform these measurements 
separately for B° and B°. A study of B° — > K+tt^tt" 
permits a measurement of the magnitude of -R3/2 but not 
its phase. Eq. I|18fl implies that I-R3/2I — 1 is proportional 
to Im(r3/2) and vanishes if r 3 / 2 is real. 

The study of B° — > Kstt + tt~, which is not flavor spe- 
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cific, is more challenging. In order to measure the relative 
phase between A(K* + n~) and A(K*~ir + ), as required by 
Eqs. J3J and 0, these amplitudes must interfere through 
B°-B° mixing leading to a common Kstt + tt~ state. Ob- 
serving this interference in e + e _ collisions at the Y(45) 
requires a time-dependent measurements using initially 
tagged B° or B° mesons. The recent time-integrated 
analysis by Belle [3^| assumed no direct CP asymme- 
try in B° — > K* + ir~ , summing over initial B° and B°. 
We note that, in fact, an untagged amplitude analysis 
does not have to make this assumption, permitting sepa- 
rate measurements for the magnitudes of A(K* + tt~) and 
A(K*~ir + ). However, measuring the relative phase be- 
tween these amplitudes requires a time-dependent mea- 
surement. 

A fourth process in this class, B + — ► Kstt + tt°, which 
has not yet been measured, determines the magnitudes of 
the four amplitudes, A(K*°ir+), A(K*+ir Q ), A(A *°7r+), 
A(Kq + 7t°), and their relative phases. Finally, a very dif- 
ficult mode which is not needed is B — > KsTT a ir a , where 
measuring the phase difference between A(K* q tt°) and 
A(K*°ir°) would require time-dependence. 

In order to apply Eq. l|2()l) . the linear constraint in 
the (p, fj) plane, where r 3 / 2 is given in (|30[) . it is suffi- 
cient to perform amplitude analyses for merely two pro- 
cesses involving neutral B decays, B° — ► K + tt~tt and 
B° — > Kstt + it~ ■ Time-dependence in the second process 
is crucial. The first process measures the magnitudes 
of A(K* + n~) and A(K* n°), their relative phase, and 
the corresponding CP-conjugate quantities, but not the 
phase difference between B° and B° decays. [Here and 
below K* denotes both if* (892) and if,* (1430)]. The sec- 
ond process measures the magnitude of A(K* + n~) and 
its CP-conjugate, and the relative phase between these 
two amplitudes. This set of measurements determining 
the complex ratio R3/2 defined in Eqs. @ and is over- 
complete since |A(A* + 7r~)| and its CP-conjugate are 
measured both in B° — » A + tt~it and in B° — > A,j7r + 7r~ . 

Charged B decays, B ± -► K ± n ± n^ and B ± -> 
Ksn^ir provide further constraints on CKM parame- 
ters using the measurable ratio Rx of B — > K*w and 
B — ► K*w decay amplitudes. This, together with the 
constraint from R3/2, leads to a highly constraining set 
of measurements for p and fj. Since the four physical 
B — > K*n amplitudes are not mutually independent [see 
the quadrangle relation Eq. (J2J], we propose studying 
B — > K*ir amplitudes in terms of the isospin amplitudes 
•B1/2, A 3 / 2 and Ax, where X corresponds to the state de- 
fined in 1|33|) . In order to demonstrate the extent to which 
these CKM constraints are over-deterministic, thereby 
permitting a precise constraint on the point (p,fj), we 
now count the number of parameters and observables. 

We have a total of eight parameters, the magnitudes 
of B1/2 and its CP conjugate B1/2, the magnitudes of 
A3/2 and Ax, the three relative phases among these four 
amplitudes, and a CKM ratio fj/(p + C). [The CP con- 
jugates ^3/2 an d Ax are not independent parameters 
and are given by Eqs. (J5J and (|26[1 .] These eight 



parameters can be used to fit seventeen observables con- 
sisting of |^(if*°7r ± | obtained from B + — > K + n + n~, 
magnitudes of ^4(A*°7r ± ) and A(K* ± tt°) and their rel- 
ative phases obtained from B ± — > Kstt^it , magni- 
tudes of A{K* + tt~),A{K* q t: ), their CP conjugates and 
their relative phases obtained from B° — > K ± tt t tt°, and 
magnitudes and relative phase for A(B° — > K* + ir~) 
and A(B° — » A*~7r + ) obtained from time dependent 
5° — > Ks7t + tt~- We have not included in this count- 
ing the decay B° — > A,g7r 7r which is most challenging. 



B. B s -» Kit it 

The weak phase 7 can be determined using Dalitz 
plot analyses for B s — > K ± tt t tt and i? s — > 
Ksir tt^ ■ These studies permit extracting the magni- 
tudes y4 s (A*+7r _ ), A s (K*°tt ), their CP conjugates and 
relative phases between these amplitudes. This leads 
through Eqs. Q-© to a measurement of the phase ^%,/ 2 , 
which gives 7 with high theoretical precision, as has been 
discussed in Section III.B. 

In contrast to the case of B° — > A's7r + 7r~ produced 
at the T(4S'), the above measurements can be performed 
with B s — > Kstt + tt~ produced at hadron colliders with- 
out the need for flavor tagging and time-dependence. 
Because of the lack of quantum coherence between B s 
and B s produced in pairs, the charge-averaged time- 
integrated decay rate for decays into a common state 
/ = Kstt + tt~ involves an interference term proportional 
to the width difference Ar s in the B s system, for which 
one expects y s = AT S /2T S = 0.12 ± 0.05 35]. 

The untagged integrated decay distribution is given by, 

d 2 T(B s -> /) d 2 T(B s -> f) _ 
dsi 2 dsi 3 dsi 2 dsi 3 

1 \(\Af + \Af)-2y s Ke(UA*)] , 

where A = A(B S -*•/), A = A(B S -> /), q/p ~ 1. 
Assuming that a reasonably accurate measurement for 
y s will exist by the time an amplitude analysis will be 
performed for this decay, the relative phase between 
A s (K* + tt~) and A s (K*~ir + ) can be measured through 
the interference term involving y s . Otherwise, a time- 
dependent measurements of this decay will be required. 

In order to show that the above relative phase is mea- 
surable using untagged B s , consider the contributions of 
A s (A*+7r-) and A s {K*-tt+) to A and A in g5J. Us- 
ing the dependence of the Breit-Wigner functions /£»+ 
and /SJl on S12 and S13, the untagged decay distribu- 
tion (1451) provides four observables (the real part of the 
interference term provides two observables) which deter- 
mine the magnitudes of A s (K* + ir~) and A s (A* - 7r + ) and 
their relative phase. While in reality this relative phase 
may be affected by interference with other resonant or 
non-resonant terms in the amplitude, this proves that, 
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once y s is given, this phase can be measured through an 
untagged amplitude analysis of B s — > Kstt + tt~ . 



C. B s -> KKtt 

As noted above, the CKM constraints following from 
amplitude analyses of B s — > KKtt decays are less precise 
than those following from studies of B — > Kirir and B s — > 
Kirn. This is due to theoretical errors in the hadronic 
electroweak penguin parameters, rf and r'{ [Eq. 1|42|) ]. 
which are larger than in r 3 / 2 [Eq. (JSUJl], rx [Eq. (j3"6*j) ] 
and r|, 2 [see discussion in Section III.B]. 

In order to obtain a CKM constraint related to the 
phase $f , for instance, one must measure the amplitudes 
in {T2J, f or B s -> K*+R- and S s -> AT* /? , their 
charge-conjugates, and the three relative phases between 
these amplitudes. This can be achieved by amplitude 
analyses for a pair of processes belonging to this class. 
For instance, using B s — > K + K~ir° one can measure the 
magnitudes of A(B S -> K*+R-), A(B S -► K*-R+), 
their charge-conjugates and the relative phases between 
these amplitudes. A study of -B s — > K + Kstt~ permits 
measurements of the magnitudes of A(Z? S — > A*°A°), 
A(B S -> K*-K+), A(B S -► K*°K°), A(B S -> K*~K+) 
and the respective relative phases. This information suf- 
fices for fixing $f. 

Decay distributions in _B S — > KKtt involve twice as 
many relevant quasi two-body amplitudes as in i? — > 
A7T7T and -B s — ► Kirn, because i? s and -B s can decay to a 
common non-flavor K*K state. The large number of am- 
plitudes and relative phases which must be determined 
in B s — > KKtt requires at least as many observables. 
While in principle possible, this seems to pose a serious 
challenge to applying this method to B s — * KKtt decays. 



V. CONCLUSION 

We have studied in great detail a method proposed in 
Ref. 0,0 for obtaining new constraints on CKM parame- 
ters using B( s ) — > (if*7r)j = 3/2 amplitudes, extending the 
method to B -> (AT*7r) J=1/2 and to B s -> if *A (.R^iT) 
amplitudes measured in _B — > A"7r7r and _B S — > KKtt, re- 
spectively. Two judiciously chosen isospin amplitudes in 
B — > A*7r have been shown to be over-constrained by 
several £? — * Ktttt amplitude analyses, providing a pre- 
cise linear constraint between the CKM parameters p and 
fj. The slope of the linear relation is a measurable quan- 
tity, while the intercept po where fj — is a calculable 
quantity involving a theoretical error of 0.03. A study of 
B s — » A*7r amplitudes in B s — > Ktttt leads to a very ac- 
curate extraction of the weak phase 7 with a theoretical 
uncertainty below one degree. 

The resulting theoretical precision in determining 
CKM parameters in B — > A7T7T and B s — ► Ktttt has been 
shown to be essentially at the level of isospin breaking 



corrections since the method is based on isospin sym- 
metry considerations, while flavor SU(3) has been used 
to estimate uncertainties from a subset of small EWP 
contributions. A larger hadronic uncertainty from EWP 
contributions is found in a CKM constraint obtained by 
studying B — > K*K and B — > K*K amplitudes con- 
tributing to B s — > KKtt. 

There is one crucial theoretical difference between ap- 
plying this method to AS = 1 B — > Ktttt and B s — > 
KKtt and applying it to AS* = B s — > Ktttt. The first 
two classes of processes are dominated by AJ = QCD 
penguin amplitudes which are eliminated in the relevant 
isospin amplitudes. In the Standard Model this implies a 
delicate cancellation between physical amplitudes defin- 
ing the numerators and denominators of the AI = 1 ob- 
servables Ri and Rf(R'f) on which the method relies. In 
contrast, in B s — > Ktttt decays the method relies on mea- 
suring the AI = 3/2 isospin amplitudes which involves 
dominant tree contributions. This would seem like a dis- 
advantage of using B — > Ktttt and B s — > KKtt relative 
to B s — > Ktttt for extracting CKM parameters. How- 
ever, this cancellation in the Standard Model turns into 
an advantage when one is searching for New Physics in 
AI = 1 operators. 

While applications of the method to B s decays can be 
foreseen in future experiments at hadron colliders, data 
for B — * Ktttt are already available from e + e~ collisions 
at the T(45*), and should be analyzed in the manner pro- 
posed here. Amplitude analyses of a few processes in the 
class B — > Ktttt have already been performed, measuring 
amplitudes and relative phases for B — > A*(892)7r and 
B -> K^(U30)tt [H IM El H|. Since the method is 
based on AI — 1 amplitudes, a first important step to- 
ward its full implementation is observing a violation of 
AI = QCD penguin dominance in these quasi two-body 
decays. 

This question has been studied recently [3(j. It was 
shown that in all cases but one AI = holds well within 
current experimental errors. For instance, AI — domi- 
nance implies 

2B{B° -» A^°7r ) = B{B° -> K*+tt-) , (46) 

which holds experimentally within large errors, in units 

of to- 6 H3, 

51.0 ± 19.8 = 46.6±|e • ( 47 ) 

The exceptional case where AJ = seems to be violated 
is the equality, 

2B{B° -» A^ 7T°) = B{B° -» K*+tt-) , (48) 

where current experimental values |22| . 

3.4 ± 1.6 = 9.8 ± 1.1 , (49) 

show a discrepancy of 3.3er. One would have to watch 
carefully whether this discrepancy holds in future mea- 
surements. 
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This method requires performing amplitude analyses 
of B — > Ktttt separately for B and B , as one must 
measure the ratio of B — ► ^*7r and £> — > .K"*7r ampli- 
tudes. The method does not require observing a direct 
CP asymmetry in Dalitz plots for B — > A'tttt or an asym- 
metry in £> — * i^*7r decay rates. We recall that no time- 
dependence is needed in order to observe direct CP vio- 
lation in the Dalitz plot of B° — > Kstt + tt~ through an 
asymmetry with respect to exchanging n + and tt~ |37|. 
We have stressed the importance of performing a time- 
dependent Dalitz plot analysis of B° — > ^571-+ 7r~, which 
is required in order to measure separately amplitudes 
for B° -> iv:* + 7r- and 5° -> i^*"7r+ and their relative 
phase. 

The method presented here for obtaining a linear re- 
lation between p and 77 in B — > _ftT7r7r may be compared 
with a study of 7 in i? — > DK |38j • The latter method 
involves an extremely small theoretical uncertainty from 
D°-D° mixing [3!j when studying CP-eigenstates and 
flavor states in D decays. Applying this method to 
non-CP and non-flavor three body D decays such as 
D a — > Kstt + it~ introduces a theoretical error in 7 caused 
by modeling the three body decay amplitude in terms of a 
sum of resonant and non-resonant contributions. Model- 
dependence in amplitude analyses for B — ► Ktttt is ex- 
pected to be larger than in D° — ► Kstt + tt~ because the 
former processes involve lower statistics and higher com- 
binatorial backgrounds. Fortunately, the uncertainty of 
modeling B — » Ktttt is mainly in non-resonant ampli- 
tudes [3(1 0^1 13 , which spread over the entire phase 
space, but less in K*tt amplitudes which are used in the 
proposed study. 

While measuring 7 from an interference of tree ampli- 
tudes in B — > DK is most likely to receive only small 
corrections from New Physics [33, E3 > the extraction of 
a linear constraint between p and 77 in B — > Ktttt may 
be affected more significantly by such effects. Thus, val- 
ues for CKM parameters obtained in the two methods 
may differ, indicating short distance b — > sqq operators 
beyond the Standard Model. The study of B — > Ktttt 



is insensitive to new AI — QCD penguin-like opera- 
tors which cancel in the ratios Ri , but is affected by new 
AI = 1 operators. Such operators are often referred to in 
the literature as anomalous electroweak (or Trojan) pen- 
guin operators . The sensitivity to such contributions 
is high because in the Standard Model AI = 1 terms in 
B — > K*tt are suppressed relative to AI = contribu- 
tions. Other tests for such AI = 1 operators have been 
proposed in terms of isospin sum rules for rates ^3 an d 
CP asymmetries in B — > Kn [43| . 

A somewhat similar situation occurs in B s decays when 
comparing the theoretically precise measurement of 7 in 
charmless B s — > Ktttt discussed here with the potent ially 
accurate measurement of this phase in B s — > D~K + |44| . 
Both methods require B s -B s mixing, but no time de- 
pendent measurement is required in B s — > Ktttt due to 
additional phase information coming from Dalitz plot in- 
terferences. In B s — Ktttt the measurement of 7 fol- 
lows from studying AI = 3/2 b — > uud tree amplitudes, 
while in B s — > D~ K + the phase occurs in the interfer- 
ence of AI = 1/2 b — > cus and b — > ucs tree amplitudes. 
Whereas New Physics operators in the latter case are pos- 
sible in principle, their effects on the determination of 7 
are less common and are expected to be much smaller 
than the effects of potentially new AI = 3/2 operators 
contributing in B s — ► Ktttt. Such AS — operators 
are usually expected in the same class of models where 
anomalous AS = 1 electroweak penguin operators occur. 
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